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Abstract 

We perform molecular dynamics simulations of a well-known water model (the TIP5P pair 
potential) and a simple liquid model (a two-scale repulsive ramp potential) to compare the regions 
of anomalous behavior in their phase diagrams. We select the parameters of the ramp potential by 
mapping it to an effective pair potential derived from the TIP5P model. We find that the regions 
of anomalous behavior in the phase diagrams of both systems can be mapped onto each other if (i) 
pressure P and temperature T are replaced by T — Tc and P — Pc, respectively, where {Tc, Pc) 
are the coordinates of the liquid-liquid critical point of the corresponding system; and (ii) a single 
ramp particle corresponds to two TIP5P molecules. We present heuristic arguments supporting 
point (ii). We also argue that the water-like anomalies in the ramp potential are due to the ability 
of the particles to reproduce, upon compression or heating, the migration of water molecules from 
the second shell to its first shell. 
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Liquid water is peculiar as reflected by its thermodynamic and dynamic anomalies 
such as the density decrease upon isobaric cooling (density anomaly) and the increase of 
diffusivity upon isothermal compression (diffusion anomaly). It has been proposed that these 
anomalies may arise from a liquid-liquid critical point (LLCP) in the deeply supercooled state 
of water Several other liquids (e.g., silica, silicon, carbon, and phosphorous) with local 
tetrahedral order ji, [s], 0] also may show water-like anomalies. These anomalies of water 
and the LLCP can be reproduced by simple liquids interacting via core-softened spherically 



symmetric potentials wnicn lacK tne strong orienti 
liquids 



lich lack the strong orientational interaction expected in tetrahedral 



Water also possesses structural anomalies which occur when metrics describing both 
translational and orientational order decrease upon compression, as found in both the SPC/E 
and TIP5P water models j^, [3]. Water's structural anomaly is also reproduced by a family 
of core-softened spherically symmetric potentials possessing two characteristic length scales 
(Jo and (Ti (see the ramp potential in Fig.[T](a)) [l5]. In order to exhibit a water-like structural 
anomaly, the ratio A = cro/o"i must lie within a small interval around 0.62, the ratio of the 
distances to water's first and second neighbor shells, 0.28 nm/0.45 nmjl^. 

A quantitative connection between the ramp potential and water's pair potential has not 
been established, as well as the relation between the regions of anomalies in their respective 
phase diagrams. In this work, we show that the effective pair potential derived from the 
TIP5P water model (20! can be approximated by a two-scale spherically symmetric repulsive 
ramp potential, allowing us to assign physical units to the temperature and pressure of the 
ramp model. We perform molecular dynamics simulations using both the TIP5P and ramp 
potentials and compare the regions of anomalies in the corresponding phase diagrams. We 
find that the regions of anomalies in both phase diagrams are quantitatively similar if (i) 
pressure P and temperature T are measured in terms of T — T^; and P — Pc, respectively, 
where (Tc, Pc) are the coordinates of the liquid-liquid critical point of the corresponding 
system; and (ii) a ramp particle corresponds to two TIP5P molecules. We present quantita- 
tive arguments supporting point (ii) and provide a simple picture to explain the similarities 
observed in the TIP5P and ramp potentials. A ramp liquid particle corresponds effectively 
to two water molecules, one molecule plus 1/4 of each of its four neighbors. The water-like 
anomalies in the ramp potential are due to the ability of the particles to reproduce, upon 
compression or heating, the migration of water molecules from the second shell to its first 
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shell. 

The TIP5P model is a well-known water model and its parameters are defined in physical 
units, so values of P and T from simulations can be compared directly with experiments 20|. 
Instead, thermodynamic properties in the ramp potential are given in terms of potential 
parameters, such as {(To,f/o}, and the particle mass, m. To compare the phase diagrams 
of the ramp potential to that of the TIP5P model, we will define ctq and Uq in units of 
'nm' and 'kcal/mol', respectively, and m in units of 'g/mol'. We do this by calculating 
Ucsi'f'), the effective spherically symmetric pair potential between water molecules from 
the T1P5P model simulations. Ucs{r) is obtained from the oxygen-oxygen pair correlation 
function g{r), by solving the Ornstein-Zernike equation and using the hypernetted chain 
approximation [21]. The resulting Ues{r) depends on T and density p |22|, but has no 
significant change for different state points in the anomalous region. For the TIP5P model, 
the range of anomalies is approximately 220 K < T < 320 K and 0.90 g/cm^ < p < 
1.16 g/cm^ [19]. We select a state point located in the middle of the anomalous regions, 
at T = 280 K and p = 1.00 g/cm^, and calculate g{r) and Ues{r) [see Fig. [l](b) and 
Fig. [H^c)]. We find that Ucs{r) is similar to the effective pair potential obtained from the 
experimental g{r) 2l|, and shows a hard- core- like steep repulsion at r ^ 0.26 nm and an 
approximately linear repulsive region covering the distance spanned by the second shell of 
a central water molecule, approximately 0.32 nm < r < 0.45 nm. The shallow minimum 
at r = 0.28 nm is caused by hydrogen-bonding attraction and corresponds to the first 
peak of g{r), while the minimum at r = 0.45 nm [Uei = t/cfT(0.45 nm)= —0.45 kcal/mol] 
corresponds to the second peak of g{r). Ucs{r) also shows a maximum at r ^ 0.32 nm 
[Ueo = f4fr(0.32 nm)= 0.66 kcal/mol] that corresponds to the first minimum of g{r). 

Figure [U^c) also shows that a ramp potential is a good approximation to Ucs{r). In the 
figure we set ai = 0.45 nm and define the ramp part of the potential such that it intersects the 
plot of Ucs{r) at {Ueo + Uei)/^. The intersection of the ramp part of the potential with the 
hard core of Ues{r) is used to define Uq and ctq. This results in o"o = 0.267 nm, which is located 
between 0.28 nm, the first peak position of g{r) and 0.26 nm, roughly the infinite repulsion 
part of Ucs{r). Therefore, A = ao/cri = 0.593 and Uq = ?7cfT(o"o) ~ f^cfr(o"i) = 1-31 kcal/mol. 
Uq is approximately the energy barrier that water molecules need to overcome to migrate 
from the second shell to the first shell positions in terms of the effective potential. It is also 
roughly the energy that ramp particles need to overcome to reach the hard core distance. 
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To define m in physical units, we argue that the ramp particle corresponds to two water 
molecules. This is based on the crystalline phases of water (hexagonal ice) and ramp po- 
tential (hep) [see Fig. [Tj . The hexagonal ice can be formed by combining units such as that 
shown in Fig. [2J Each of this units is formed by a central water molecule and 1/4 of its four 
neighbors, which are tetrahedrally arranged. To form the hexagonal ice, such units must 
form an hep network. Therefore, if the crystalline structure of the ramp potential model is 
identified with that of hexagonal ice, a ramp particle must be identified, on average, to the 
unit shown in Fig. [2l The mass of a water molecule is = 18 g/mol, thus, the mass of 
a ramp particle is m ^ (1 + 4 x 1/4) m^; = 36 g/mol. Alternatively, the present argument 
implies that the number density of the ramp potential model corresponds to twice the num- 
ber density of water, and this will be relevant when comparing the pressures of the ramp 
and TIP5P models [23|. To test the idea that a ramp particle corresponds approximately 
to two water molecules, we calculate the average number of neighbors, Nq, that a water 
molecule has within a distance of r < (Tq = 0.267 nm. Using the g{r) from Fig. [T] we get 
A^o = Airn r''^g{r')dr' ^ 1 (here, n is the number density), in agreement with our view. 
The correspondence between one ramp particle and two water molecules is also supporte d by 



computer simulations of the ramp potential with A = 0.581 and an attractive part 
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li. 



Such a ramp potential model has both liquid-gas (LG) and liquid-liquid (LL) critical points. 
Application of the values for ctq, Uq, and m that we use here to the data from [l4| results 
in plg ~ 0.314 and pn ^ 1.188 g/cm^. These values approximately coincide with the 
experimental critical den sity of water pic ~ 0.322 2^ and the LL critical density of TIP5P 
water pll ~ 1-13 g/cm^ |25 |. 

To compare the regions of anomalies in the phase diagrams of the TIP5P and ramp 
potentials, we first obtain the LLCP coordinates, {Pc,Tc,pc)- The LLCP in the TIP5P 
model is accessible in MD simulations and is located at Tc = 217 K, Pc = 340 MPa, and 
Pc = 1.13 ± 0.04 g/cm^ [261, l27|]. Instead, for the ramp potential of Fig. [I](a), the LLCP is 
located at temperatures below those accessible in simulations 13]. In this case, the LLCP 
can be located by extrapolating the isochores in the P — T phase diagram to low-T (the 
isochores cross each other at the LLCP). This procedure indicates that the LLCP is located 
at Tc = 16.5 K, Pc = 967 MPa, and p^ = 1.19 g/cm^. 

Figure [3] shows the phase diagrams of the TIP5P and ramp potential models, obtained 



by MD simulations (for details see |l2l. Il5l. Il9|). To emphasize the quantitative similarities 
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of these diagrams we place the origins of P and T axes at the LLCP of the corresponding 
models. In both models, the density anomaly region is within the diffusion anomaly region, 
which is enclosed by the structure anomaly region. A comparison of panels (a) and (b), or 
(c) and (d), shows quantitative similarities in the regions of anomalies of both models. The 
density anomaly region covers approximately the ranges —500 < P — Pq < MPa, T — Tc< 
60 K, and 0.95 < p< 1.15 g/cm^, for the TIP5P model, and -520 < P - Pc < -10 MPa, 
T — Tc < 70 K, and 0.85 < p < 1.20 g/cm^ for the ramp potential model. Similarly, 
the diffusion anomaly region covers approximately the ranges —500 < P — Pc < 50 MPa, 
T - Tc < 90 K, and 0.92 < p < 1.17 g/cm^ for the TIP5P model, and -520 < P - Pc < 
MPa, T — Tc < 90 K, and 0.85 < p < 1.21 g/cw? for the ramp potential model. The 
structure anomaly region, defined by the loci of the order parameter extrema (tmm, Imax) and 
{imin, Q&max) for the TIP5P and ramp potential models respectively, covers approximately 
the range -500 < P - Pc < 150 MPa and 0.9 < p < 1.22 g/cm^ for the TIP5P model, 
and —500 < P — Pc < 150 MPa and 0.85 < p < 1.27 g/cm^, for the ramp potential model. 
However, the shape of the structural anomaly region is different in these models. It expands 
up to T — Tc < 100 K for the TIP5P model, while for the ramp potential it expands to 
much higher temperature outside the graph. 

A possible reason for the quantitative similarities in the regions of anomalies of water 
and ramp potential model is that this model is able to reproduce quantitatively the observed 
migration of water molecules from the second shell toward the first shell upon compression 



or heating [19|, |28|, |29|]. We discuss first the probability distribution, P(ri), of the distance 



between a ramp particle and its nearest-neighbor. Figure Hl^a) shows the evolution of P(ri) 
upon isothermal compression. As density increases, the maxima of P(ri) shifts from r = 
0.42 ^ (Ti, at low density, to r = 0.267 = ctq, at high density. Figure|l](b) shows that a similar 
but less pronounced changes in P(ri) occur upon isobaric heating. Thus, upon compression 
or heating, particles move from the soft-core distance (corresponding to water's second shell) 
toward the hard-core distance (corresponding to water's first shell) of the ramp potential. 



Similar structural changes occur in water [19|, l28|, l29|] . In particular. Figures HTa) andlD^b) 



can be compared with the corresponding Figs. 4(c) and 4(f) of ref. 19|] obtained for the 
TIP5P model. The probability distribution, P{Qq), of the orientational order parameter, 
Qe 13 5 of the ramp potential particles is shown in Figs. |l](c) (upon isothermal compression) 
andlD^d) (upon isobaric heating). Upon compression or heating, the maximum P{Qq) shifts 
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to small values of Qe, i.e. orientational order decreases. Similar structural changes occur 
in the orientational order of water's second shell; Figs. Hl^c) and IH^d) can be compared 
with the corresponding Figs. 4(b) and 4(e) of ref. 19|] obtained for the TIP5P model. For 
a quantitative comparison of the structural changes in the ramp and TIP5P models, we 
calculate the number of neighbors, N{r), as a function of the distance r from a central 
water molecule/ramp particle in both models. The increase of N{r) with density, AN{r), is 
shown in Fig. Hl^e). We see that AN{r) for both models are remarkably similar and overlap 
for approximately r > 0.37 nm. Thus, the ramp potential reproduces quantitatively the 
migration of water molecules from the second shell toward the first shell upon compression 
or heating. 

In summary, our study makes a microscopic quantitative connection between a ramp po- 
tential and TIP5P water model and shows that orientational interactions, such as hydrogen 
bonding, are not necessary to reproduce water-like anomalous properties. In general, the 
ramp potential provides an understanding of the anomalous features of tetrahedral liquids. 
These features are caused by a large empty space around the tetrahedrally coordinated 
molecules, which is reduced as temperature and pressure increase. In the ramp liquid, this 
empty space is created by the repulsive soft core. 
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FIG. 1: (Color online) (a) The two-scale spherically symmetric repulsive ramp potential: do corre- 
sponds to the hard core distance, ai characterizes a softer repulsion range that can be overcome at 
high P and T. The central ramp particle (black) and its twelve nearest neighbors (yellow) form a 
hep crystal structure in a range of densities corresponding to the density anomaly. In our study the 
parameters are defined in physical units as ai = 0.45 nm, ao = 0.267 nm and Uq = 1.31 kcal/mol 
by mapping the ramp potential to the effective potential in (c) [see text], (b) The pair correlation 
function, g{r), and (c) spherically symmetric effective potential, Ucs{r), from the simulations using 
the TIP5P model at T = 280K and p = 1.00 g/cm^ (solid line). For hexagonal ice, the twelve 
neighbors (yellow) in the second shell of the center water molecule (black) also has a hep structure 
while the four nearest neighbors (green) in the first shell are located in the corner of a tetrahedron. 
Ucs{r) can be approximated by a ramp potential (dashed lines). By calculating the area of g{r) for 
r < o"o we find that the hard core of the ramp particle roughly incorporates two water molecules 
(see also Fig. [2]). 
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FIG. 2: Sketch of a water molecule and 1/4 of each of its four nearest neighbors in a tetrahedral 
arrangement. Only oxygen atoms are shown for clarity. Each of these units corresponds effectively 
to a ramp particle. The hexagonal ice (the low pressure crystal of water) can be obtained by 
combining these units in an hep lattice (the low pressure crystal of the ramp potential model). 



We notice this figure is in the spirit of Walrafen pentamer 
consists of only two water molecules. 
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FIG. 3: (a) Three anomalous regions of TIP5P water in modified P-T phase diagram. The dashed 
lines are the isochores with density />=1.20, 1.16, 1.12, 1.08, 1.04, 1.00, 0.96, 0.92, 0.88 g/cm^ 
from top to bottom. Density anomaly region is defined by TMD (temperature of maxima density) 
lines, inside which the density increases when the system is heated at constant pressure. Diffusion 
anomaly region is defined by the loci of DM (diffusion maxima- minima) , inside which the diffu- 
sivity increases with density at constant T. Structural anomaly region is defined by the loci of 
translational order minima (tmin) and maxima (tmax)) or orientational order maxima gmax {Q&max 
for ramp liquid), inside which both translational and orientational orders decrease with density 
at constant T (see refs. 15, [2] for details). Here t quantifies the tendency of molecular pairs 
to adopt preferential separations, and q quantify the local tetrahedrality of water (Qg quantify 
the local orientational order of twelve nearest neighbors in the first shell of a ramp particle), (b) 
Anomalous regions for the ramp liquid, here the values of P and p are doubled in order to compare 
with the corresponding values of water [see text]. The dashed lines are the isochores with density 
p=1.33, 1.28, 1.23, 1.18, 1.14, 1.09, 1.05, 1.02, 0.98, 0.94, 0.91, 0.88, 0.85, 0.82, 0.79, 0.77, 0.74 
g/cm^ from top to bottom. The structural anomaly region is open and converge at much higher T 
compared with water, (c) and (d) are the anoilig,lous regions in the T-p phase diagrams. 
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FIG. 4: Probability distribution of the distance ri of a central ramp particle and its nearest 
neighbor at (a) constant T and (b) constant p. (c)-(d) Probability distribution of orientational 
order parameter for a ramp potential particle corresponding to panels (a) and (b), respectively. 
Upon heating or compression, ramp particles move from the soft-core distance toward the hard- 
core distance and the orientation order parameter decreases. Similar structural changes occur in 
water [l^ . (e) Increase in the number of neighbors, AA^(r) = N{r)\p^ — A^(r)|pp, where pi = 
0.88 g/cw? and po = 1-08 g/cm'^, for the TIP5P and ramp potentials. AA^(r) in both models are 
remarkably similar and overlap for approximately r > 0.37 nm. Since a ramp particle corresponds 
to two water molecules, we doubled the values of N{r) and AA^(r) obtained from the simulations 
using the ramp potential model. 
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